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Abstract 

We study the corner-to-corner resistance of an M x iV resistor net- 
work with resistors r and s in the two spatial directions, and obtain an 
asymptotic expansion of its exact expression for large M and N . For 
M = N,r — s = l, our result is 

4, 0.266 070 0.534 779 ^,1, 
Rm.n = - logiV + 0.077 318 -f —^^ + 0{—^). 



1 Introduction 

A classic problem in the theory of electric circuits is the computation of 
the resistance between two nodes in a resistor network. Formulated by 
Kirchhoff [1] more than 160 years ago, the problem has been studied by 
numerous authors over many years (see, for example, [2l[3]). Kirchhoff 
explored the graph-theoretical aspect of the algebraic formulation and 
obtained the two-point resistance in terms of 2-rooted spanning forests 
and spanning trees. But the formulation, while elegant, does not pro- 
vide sufficient physical insights. Past studies have instead focused on 
infinite networks for which analysis can be carried to fruition [3]. 

The computation of the asymptotic expansion of the corner-to- 
corner resistance of a rectangular resistor network has been of interest 
for some time, as its value provides a lower bound to the resistance 



of compact percolation clusters in the Domany-Kinzel model of a di- 
rected percolation [5]. The corner-to-corner resistance has been stud- 
ied by one of us (JWE) numerically using the method of a differential 
approximants [5] together with a Neville table analysis [7]- 

Recently, one of us (FYW) has re-visited the two-point resistance 
problem [8], and deduced a closed-form expression for the resistance 
between arbitrary two nodes for finite networks. However, the exact 
expression obtained in [S] is in the form of a double summation whose 
mathematical and physical contents are not immediately apparent. In 
this paper, we take a closer look at this summation formula and obtain 
its asymptotic expansion for large lattices. 

The organization of this paper is as follows: In Sec. 2 we recall 
the expression of the corner-to-corner resistance in an M x N resistor 
network obtained in [8], and reduce it to a form more manageable for 
our purposes. One of the two summations in the resistance expression 
is carried out in Sec. 3 by using a new summation identity which we 
derive. The resulting expression is written in the form of a dominant 
term plus a correction. Asymptotic expansions of the dominant and 
correction terms are obtained in Sees. 4 and 5, and we summarise the 
results in Sec. 6. We also show that the exact expression of the asymp- 
totic expansion is in agreement with those determined numerically [7] . 



2 Formulation of the summation formula 

Consider a rectangular M x N network of resistors with resistances r 
and s on edges of the network in the respective horizontal and vertical 
directions. For definiteness, we consider both M, N even, and expect 
the asymptotic expansion to be independent of this choice. The exam- 
ple of an M = 6, iV = 4 network is shown in Fig. 1. 

(M-1,N-1) 



(0,0) r r r r r 



Figure 1: An M x N resistor network. 



Using Eq. (37) of [S] , the resistance between opposite corner nodes 
(0, 0) and (Af - 1, TV - 1) of the network is 

/ N r(Af-l) .s(A^-l) 
R{MxN}{r, s) = + — 

2 ^i^i ^1 [cos cos (i0„) - cos (M - cos {N - i),/), 

^ — E 



MN ^ ^ r^ifl - cos6',„) s-i(l - cos0„) 

m— 1 n— 1 ^ ' ^ ' 



(1) 



2 



where 9m = Tmr/M,(f)„ = m:/N. Re-arranging the numerator in the 
summand, H]) becomes 



RMxN{r,s) 



r{M - 1) s{N - 1) 



N 



Af-l N-l 

EE 



M 

cos2(0„/2)cos2(0„/2) 



MN ^ r-i(l-cos6l„)+s-i(l-cos(/.„) 

m=l n=l(m+nodd) 



(2) 



There are two possibihties for the restriction m -f n = odd to hold, 
namely, 

m = 239 - 1, n = 2q, _p = 1, 2, .., M/2, (7=1,2, iV/2, 
n = 2p- l,m 2g, p = 1, 2, .., iV/2, g 1, 2, Af/2. 

Splitting the sum into two parts accordingly and introducing the no- 
tation 

A-^, ^^=1^-2Jm' 
we obtain 

RMxN{r,s) = {rs)i[RMxN{r/s) + RNxM{s/r)] (3) 

where 

- 1) 



Rmxn{p) 



AT 

A//2 Af/2 

EE 



COS^ Aq{\+ pSlT? Aq) 3 x 

9 A 9 COS -/Ifj 

p sin + sin Bp 



(4) 



p=l q=l 

Sums of the term cos^ Aq can be carried out using the identity 
}_^cos ( 

q=l 

This yields 



4 2 



(5) 



and 



Rmxn{p) = ^fp ^ ^ ) + '5'MxAr(p) 



RMxNir, s) = y/rs 



+ SMxNip) + Snxm{'^/ P) 



(6) 



where 



4 ^ 

Smxn{p) = -jf Y.^cos' Aq){l -f psin^ Aq) Sq ,aln{p) (7) 



q=l 



3 



with 



A//2 



(Af/2)-l 



k=0 



(8) 



3 Evaluation of Sq^M,N{p) 



It is tempting to evaluate the summation ^ by using the Euler- 
Maclaurin summation formula. But as shown in the Appendix the 
Euler-Maclaurin summation is inadequate since it does not determine 
an error term which cannot be ignored. We proceed here to evaluate 
Sq,M,N{p) by using a summation identity which we state as a lemma: 



Lemma: 

{M/2)-l 

V - 



-R{v* 



Mtanh(7r?/*) 



fc=0 

2 y/p sin Aq^l + p sin^ Aq 

where M = even and y* — y* jv(/°) defined by 

ny* 



(9) 



sinh ■ 



M 



=^sin Aq. 



(10) 



x = 



_ M 




y = 



Figure 2: Contour of integration C in (jlip . Solid circles denote simple poles 
enclosed by C . 



Proof. Consider the contour integral 



, , 1 / 7rtan(7rz)(iz , , 

Jq,M,N{p) ^ TT- f . 2 rnz^ ■2a H) 

2Tn Jc sm^ {^) + psm^ Aq 
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where the contour C consists of the lines 



M 



-oo, y ■ 



oo 



(12) 



and the imaginary axis x ~ with two half circles of radii e — > 
around the two points z — ±iy* as shown in Fig. 2. The contour 
encloses ^ + 2 simple poles of the integrand at z = ± iy* and z = 
^, |, ^^^Y^- The residue is R{y*) at the simple poles on the y-axis 



0,1,. 



and -[psin^Ag + sin^(A:+ i)^]"^ atz = A:+i, 

The integration along the contour C vanishes on the lines y = ±oo, 
and on the straight line portions of a; = 0, ^ since the integrand is odd 
in y. Hence the contour integral is nonzero only on the two half circles. 
The integrand is odd in z so that the integral along the lower half 
circle is equal to the integral in the anti-clockwise direction along the 
reflection of the upper half circle in the y— axis. The integral Jq^M,N{p) 
along the contour C may therefore be obtained by integrating round a 
circle centered on iy* . Thus, by the residue theorem, the residue at iy* 
is equal to the sum of the residues of the + 2 simple poles enclosed 
by C, hence 



(Af/2)-l 



Riy*) = 2R{y*) - V — ^ . ^ , 



which yields ©. 



(13) 
□ 



The substitution of © into ^ and ^ now yields 

tanh(7r?/*) 



^q,M,N 



(P) 



2 y/p sin AqJ \ + p sin Aq 



(14) 



and 



where 



N/2 

Smxn{p) = Dq,N{p) tanh(7ry*), 

9=1 



Dq,N{p) 



2 cos^ AqJl + p sin Aq 



N 



sin A„ 



(15) 



(16) 



Anticipating that the dominate contribution of Smxn{p) is given 
by with tanh(7r?/*) replaced by 1 (see Appendix), we rewrite 
as 



where 



SMxNip) = S'}^\p) + Am,Jv(/3), 



N/2 

4^)(p)=^i^,,„(p) 
9=1 



(17) 



(18) 
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is the dominate contribution, and 

^m,n{p) = \m,n(.p) (19) 

9=1 

is the correction with 

\m.n{p) = I?,^Ar(p)[tanh(7ry*) - l] . (20) 

Numerical evaluation of the difference Ag_M,A'(l) using tanh(7ry*) 
given by (O for M = and small values of q shows that it initially 
decreases with N but ultimately shows a rapid increase. For q = 1 the 
turning point is iV = 6 and for g = 2 it is = 12. However Aq.M,Jv(p) 
for fixed N decreases exponentially with increasing q, a fact which will 
be seen to hold for general M and N later (see Eq. ([55]) below). The 
sum in therefore converges rapidly. 

The two terms s'}^\p) and Am, Nip) in (HZ]) are evaluated in the 
next two sections. 



4 Evaluation of S'}^\p) 

The asymptotic form of S^\p) given by the summation (fT5|l is now 
deduced using the Euler-Maclaurin sum formula ([3] equation 5.8.13) 



E /p = T r f{^)dx + \ [f{Xr) - /(Xo)] 
p=l •'^0 ^ 



+ E ^tI^ [/''^-^H-.) - + EM (21) 

where /p is such that fp = /(xq +p/i), the integer r is finite and the 
error term is given by 

n L2m+2 

Erain^) - r /^^'"+^)(77^), xo < ry™ < x.. (22) 

(2m + 2)! 

But the direct application of (PT|) to effect the summation in leads 
to a divergent integral so we add and subtract 1/Aq to the summand 
and use ([2T|) with /(x) given by 



/(a:) ^ /p(x) = £^7l + psin^ X - i. (23) 
sm X ^ X 

Using xq = Q,Xr = 7r/2, /i — T:/N,r = N/2 and since fp{x) does not 
diverge at small x, the error term Em is of the order of 0{N~'^'^™^^) 
and can be neglected in m ^ oo. Denoting by Un{p) and Lm{p) the 
respective correction to the integral at the upper an lower limits, we 
obtain 

S^n\p) = I{P) + Sn + Un{p) + Ln{p) , (24) 
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where I{p) is the integral 



Hp) 



ir/2 



fp{x)dx 

1 + 4 log 2 - 2 log TT - log( 1 + p) 



p-1 
VP 



tan ^ 



(25) 



The second term in (j24l) is the added summation S'jv, which can be 
evaluated using the result ([5] chapter 5, problem 26) as 



, 2^^ 1 2^^1 2. N 1 ^ 

q=l q=l ^ m=l 



2m 



2miV2m 



(26) 



where 7 = 0.577 215 664 901 53 . . . is Euler's constant. 

The first part of fp{x) is antisymmetric about 7r/2 so the odd deriva- 
tives at the upper limit arise entirely from the — l/x term and is inde- 
pendent of p. Hence for j odd f^/\-K/2) = {-iy+^ j\{2/ Tiy+^ and the 
correction to the integral from the upper limit is 

.2^-l _2 2 



N 



2 \ ^ -822^ 



N 



4=1 



jp 



nN ' TT 2iN^i 

i—1 



E 



(27) 



which, as m — *■ cxD, cancels terms of the inverse powers of N in Sn- 
At the lower limit we have /p(0) = and 



(28) 



Using Bernoulli numbers B2 = 1/6,^4 = -1/30,56 = 1/42 (0 equa- 
tion 5.8.8), the leading terms in Ln are 



LNip) 



12iV2 



720A^4 



/f (0) 



- f(5)(0) + O(— )" 

302407V6-''' ^""^^"^^^s^J 

(29) 



with 



4'H0) = ^(-5 + 3p), 

(0) = ^(67-210p-45p2), 



1 



/^(O) = Y^(-95 + 3843p+ 2835p^ + 945p^). 
Combining ([24|) - ([27]) . we obtain the result 



(30) 



N 



(31) 
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5 Evaluation of Aq^M,N{p) 

We now evaluate A m,n (fi) given by the summation p9|) with Aq_M,JV {p) 
given by PP)) . 

For large M, with M/N — A fixed, we use 



sinh "'^(y^sinx) = y/px 



1 + P 2 , (l + p)(l + 9p) 



6 



120 



and (fTO|) to obtain 

^ ^ ^-^L 6 viv/ 120 vTvy 

where q = X^fpq. This leads to 
tanh(7rj/*) = tanh(7rg) ^-^(tt (j) sech^(7r ^-^^ 



(32) 



(33) 



^(l+p)(l + 9p) 
36 



■(1 + p)^tanh(7rg) 
Substituting into (PU]) . we obtain 
^q,M,N{p) = L'q,7v(/5)[tanh(7rg) - l] 



sech' 



(34) 



q,N 



(P) X (TTf?)!- 



+ (1 + p)sech^(7rq) 
Rewrite Dq,N{p) given by (fT6|) as 



l + 9p + ^ (7rg) tanh(7r(7) 



120 



36 



(35) 



-[2 + 2/W(0) 



^/f(0)(f)V---]. (36) 



where the derivatives are given in pop . This leads to the desired 
asymptotic expansion 



Aq,M,jv(p) = X! 



Ag,2i(A, jO) 



(37) 



with expansion coefficients 
A5,o(A, p) = — [tanh(7r(7) - l] , 



A,,2(A,p) = 2^q/W(0)[tanh(7rg) - l] - + p)sec\,\^q), 

(0) 



Ag,4(A,p) 



■ [tanh(7r(7) — l] 



\^p[TTqf[l + p) 



53 - 3p (1 + p) 



180 



18 



(38) 

(ttq) tanh(7rq) sech^(7rq). 
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As remarked earlier, values of these coefBcients decrease exponentially 
as q increases. 



6 Results 

6.1 Summary of asymptotic expansions 

Results obtained so far may be summarised as follows: the resistance 
RMxN{r, s) is given by ([3]), with Rmxn{p) expanded as 

Rmxn{p) - - logiV + - + - [27 - 1 + 2 log (- 

_ ^ N/2 
- log(l + p) + ^— tan^i ,/^] + Ln{p) + J2 \.m,n{p), (39) 

^ q=l 

where 7 = 0.577215 664 90153. . . is Euler's constant, Li^{p) is given 
by dSH]) and A,,m,w(p) given by jST]). 

As TV CX3 with A = M/N fixed, ([39|) can be written as 

RM.Nip) = ^ logiV + C(A, p) + ^Hi^ (40) 



where 



C(A, p) = Vp(a - i) + i [27 - 1 + 2 log (^) 

- log(l + p) + tan-i x/pj + ^ Ag,o(A, p), 



q=l 



&2.(A,p) = -( ^^g'p )/f-^)(0) + ^A,.,.(A,p). (41) 



(2z) 



9=1 



Here, the Bernoulli numbers are B2 — 1/6,-64 = —1/30,-66 — 1/42 
([9] equation 5.8.8). The function fp{x) is defined by ((23|) and its first 
few derivatives are given in (|30p . Equation (|37p gives an expansion 
of Aq^M,N{p) in inverse powers of A^-^ correct to 0{1/N'^) and the 
coefficients decay exponentially with q so that accurate results may be 
obtained using only the first few terms of the sum. This is illustrated 
in Table [T] in the case A = p = 1. 

6.2 The case M=N, r=s=l 

For an A^ X A^ network with ?- = .s = 1 we have X = p = 1. From ([3]) 
and (PP]) we obtain 

Rnxn{^A) = 2-RArxAf(l) 

= ^°g^ + ^" + ^ + |r + o(]^)' (42) 



9 



<? 


2Aq,o 


2A,,2 


2A,,4 


1 

2 
3 
4 
5 


-0.0047465399754997281 
-4.4402067094342628 10"*^ 
-5.5279070728467383 10"'' 
-7.7422874638854272 IQ-^^ 
-1.1566622761121781 10"" 


-0.082316647898659221 
-0.00067582947581056974 
-2.9215219029290850 10"'' 
-9.8350012986547643 10"'' 
-2.8935175541424704 10"" 


0.038515173969807909 
-0.032940604383097552 
-0.00060979439982744305 
-5.3264158004237130 10"^ 
-3.2098297733739912 10"** 




-0.0047509857178701073 


-0.082995408760387631 


0.004959416517708477 



Table 1: The coefficients Ag,2i(l, 1) in (pj) . 



where 



CO = 2C(l,l) + 2^Ag,o(l,l) 



9=1 



1 + ^ 
TT 



4 ^-^ /tanh(7r(7) — 1 



9=1 



= (0.082 069 879 627 328 • • • ) - (0.004 750 985 717 870 046 5 • • • ) 
= 0.077 318 893 909 458 - •• , 

OO 

C2 =-27rB2/l'H0) + 2^A,,2(l,l) 

9=1 

= 0.266 070 441 638 478 • • • , 



C4 



^/f)0)-f2^A,4(l,l) 

9=1 



6 



= -0.534 779 473 843 066 • • • 



(43) 



where we have used the data in Table [T] This reproduces numerical 
values of the coefficient cq determined from a differential approximant 
analysis [5] of the first 29 values of RpfxN{^, 1) together with a Neville 
table analysis [7]. Note that the correction to the dominant contri- 
bution in Co is not negligible. We have further extended the Neville 
table analysis of [7j to the next two coefficients, and obtained results 
in agreement with the theoretical values of C2 and C4. 

Finally, the asymptotic expansion (|42p is to be compared to that 
of the resistance between nodes (0, 0) and {N —1,N—1) in an infinite 
square lattice [4], 



Rn X N, oo{^, 1) = - log + 7 -I- 2 log 2 

TT 



(44) 
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Appendix 



In this Appendix we evaluate Sq^M,N{p) given by the summation 
using the Euler-Maclaurin sum formula ([S] equation 5.8.18) 

X^ff/c+i = T g{x)dx 



(45) 

where 5(0;) is such that gi = g{x{) + i/i), the integer r is finite, and 

n _9^1^2m\r> L2m+2 

(Zm, + 2)1 

where B2m are Bernoulli numbers. 

Using gS]) with xo = 0, Xr = 7r/2, h = tt/M, r = M/2, 

9i^) = .2/^.2 . (46) 
p sm + sm X 

and noting that the odd derivatives vanish at the endpoints, we obtain 

1 r^' 

Sq,M,N{p) = - g{x)dx + E„i{q, M, N, ^,„) 
I" Jo 

^ +i;™(q,M,A^,^™). (47) 



2^ sin y 1 + p sin Aq 

Comparison of (|T7|) with indicates that the dominant contribution 
of Sq^M,N{p) is precisely ([14]) with tanh(7ry*) replaced by 1, a result 
we quoted earlier. It also identifies the error term to be 

M, N, U = '--^(-f)~\ , (48) 

2y/p sin W 1 + p sin^ Aq 

a result which cannot be deduced from the Euler-Maclaurin formula. 
We point out that since the denominator of (|46p can be very small for 
q and x small, Em{q, M, N, ^,„) does not necessarily vanish even in the 
limit of m — > 00 . 
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